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Abstract 

We obtain a novel result in QCD for long range rapidity correlations 
between gluons produced in the collision of saturated high energy hadrons 
or nuclei. This result, obtained in a high energy factorization framework, 
provides strong justification for the Glasma flux tube picture of coherent 
strong color fields. Our formalism can be applied to "near side ridge" 
events at RHIC and in future studies of long range rapidity correlations 
at the LHC. 



1 Introduction 

Long range rapidity correlations in high energy hadronic collisions are of inter- 
est in QCD because causality dictates that these correlations are produced at 
very early times. They therefore provide insight into how color correlations in 
the hadron wavefunctions become dynamically manifest in multiparticle final 
states. Recent observations in nucleus-nucleus collisions at RHIC of a "near 
side ridge" structure in two-particle correlations [1|2|3| and significant forward- 
backward multiplicity correlations have rcinvigorated interest in the underly- 
ing dynamics of these correlations in QCD. At the LHC, with its wider rapidity 
coverage, such correlation studies can prove a powerful diagnostic tool both 
of multiparton correlations in QCD and of highly coherent strong color fields 
generated at early times in nuclear collisions. 

Long range rapidity correlations were previously studied in color flux tube 
models where the non-perturbative dynamics is at the QCD scale A^^^ ~ 1 
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fm~^. However, the rapid growth of parton distributions, and the requirement 
that occupation numbers in QCD saturate at l/ccg, where as is the QCD 
coupUng constant, suggests that the dynamics of color correlations is controlled 
instead by a semi-hard saturation scale Qs 3> ^qcv The properties of glu- 
ons with maximal occupation are described by the Color Glass Condensate 
(CGC) [S]; this saturated regime and the approach to it can be computed in a 
weak coupling effective field theory (EFT). 

In the CGC EFT, partons with longitudinal momenta fc+ larger than a cutoff 
momentum A"*" in a high energy hadron or nucleus (moving in the +z direction) 
are described as static light-cone color sources while modes with fc+ < A+ are 
treated as QCD gauge fields that couple to these color sources [B]. Because the 
physics is independent of this separation scale, one obtains a renormalization 
group (RG) equation-the JIMWLK equation [7] -describing the evolution of the 
distribution of fast sources as A+ is lowered. 

The QCD matter formed immediately after a nucleus- nucleus collision is 
a Glasma [51. At leading order (LO), solutions of the Yang-Mills equations 
reveal that the Glasma corresponds at early times r < to highly coher- 
ent longitudinal chromoelectric and chromomagnetic field configurations |9ll0j 
with maximal occupation numbers l/ag. At leading order, the Glasma fields 
are invariant under boosts in the space-time rapidity rj and only depend on 
their transverse positions in the nuclei and the proper time r. Further, the LO 
Glasma fields have the spatial structure of flux tubes stretching between the 
two nuclei, each localized transversely in a region of size 1/Qs- This geomet- 
rical picture naturally explains the rapidity correlations observed in the near 
side ridge in heavy ion collisions |ll|12j . However, the boost invariance of the 
Glasma configurations at LO is broken by quantum effects and it is important 
to understand their impact on multiparticle correlations. 

In two previous papers [13114] -hereafter referred to as papers I and II- we 
applied this effective field theory to the inclusive multigluon spectra in nucleus- 
nucleus collisions. The main result in these papers is a proof of the fact that 
all the leading logarithms that arise in loop corrections to these quantities can 
be absorbed into universal distributions for the fast sources of the two nuclei. 
However, in the case of the multigluon spectra, our proof was limited to the very 
peculiar situation where all the observed gluons lie in a small region in rapidity 
(of size Ay <^ as~^). This limitation was due to the fact that we did not resum 
corrections of the form as | j/i — yj \ where yi.j are the rapidities of the gluons i and 
j. These corrections become important when the rapidity separation between 
the observed gluons is of order as~^ or larger. Physically, these corrections arise 
from the radiation of extra gluons between those that are measured. This has 
a high probability of occurring if the rapidity interval between two measured 
gluons is larger than as~"^. 

The goal of the present paper, the third in this series, is to extend the treat- 
ment in the previous papers to compute inclusive multigluon spectra (to leading 
logarithmic accuracy) for arbitrary rapidity intervals between the observed glu- 
ons. In the case of the two-gluon spectrum, this is the basis for a detailed quan- 
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titative study of long range rapidity correlations in heavy ion collision^ These 
results are new and are valid (for any number of colors including the physical 
Nc = 3) in a weak coupling scheme where higher order in as contributions- 
enhanced by the same powers of the rapidity-are resummed to all orders. In 
this leading logarithmic approximation, we will demonstrate that expectation 
values of operators can be factorized as a convolution of density functionals 
from each of the nuclei times the operator computed with leading order classi- 
cal fields. These density functionals evolve according to the JIMWLK equation 
and are universal; they can be extracted, for instance, in electron-nucleus col- 
lisions. Albeit the focus here will be on nucleus-nucleus collisions, the results 
also apply to the collision of two "saturated" hadrons at very high energies. 

We note that while multiparticle correlations in the strong interactions have 
been extensively studied [15ll6j , none of the literature addresses nucleus-nucleus 
collisions, for finite Nc, in a framework where gluon fields are the dynamical 
degrees of freedom. We will comment later on interesting earlier studies [15] on 
two-particle correlations in the context of Local Reggeon Field Theory. 

The essence of the "technology" needed to resum all the leading logs in 
multigluon spectra for arbitrary rapidity separations is already contained in pa- 
pers I and II, albeit in a somewhat hidden form. In section[21 we review the main 
results of these papers and we prove a general formula for the renormalization 
group flow in the CGC when one moves the cutoff A+ of the effective theory 
by an infinitesimal amount. In section [31 we show how 1- and 2-gluon inclusive 
spectra can be obtained from this master formula. Our formula for the 2-gluon 
spectrum is expressed in terms of the usual distributions of color sources, and 
of a new object that has the interpretation of a propagator (in functional space) 
for the JIMWLK evolution. In the limit where the two gluons are nearby in 
rapidity, we recover known results for the 2-gluon spectrum. Finally, in section 
[31 we discuss the complications that arise when one wants to extend these re- 
sults to the case where one (or both) of the projectiles is dilute instead of being 
in the saturation regime. We end with a brief summary and outlook for future 
work. 

2 Leading logarithms at NLO 

Consider an inclusive multigluon field operator O. In papers I and II, we 
have shown that its LO value O^^ can be expressed in terms of light-cone 
gauge classical Yang-Mills solutions A'^, with retarded boundary conditions 
lima;0^_tx) A'^ix) — 0. In coordinate space, in the CGC effective field theory, the 
classical sources are localized along the light-cones in two strips < x~ < 1/A+ 

^In |11I12| . a simpler leading order formula was used that does not resum the corrections in 
aalyi — Ujl- This was sufficient to justify the existence of long range rapidity correlations and 
to suggest its relevance for the RHIC data. In particular, for STAR data, where the relevant 
rapidity window is Ay ~ 1.5, these resummation effects are not likely to be large. Resum- 
mations of long range rapidity corrections are however expected to improve the quantitative 
description of the STAR and PHOBOS data |2I3| and in future of the LHC data that has a 
significantly wider rapidity coverage. 
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Figure 1: Initial space-time surface S used in the initial value problem for the 
retarded classical field A^. 



and < a;+ < 1/A^. We denote by E the surface located at a distance 
Sx^ = l/A^ above the backward light-cone, as illustrated in fig.[T] Because the 
classical fields involved in obey retarded boundary conditions, O^q can be 
obtained by solving an initial value problem with initial conditions defined on 
the surface E. 

It is convenient to describe color sources in the nuclei by distributions of 
Wilson lines 

1 

^i.2{y,x±_) =T expig dz^ —^pi^2iz^ ,x±) , (1) 
JO V j_ 

where pi^2 are the color source densities in Lorenz gauge. Here the upper bound 
in the integral is related to y by y = ln(P^a;^), with the longitudinal 

momenta of the respective nuclei 0. 

In papers I and II, the next-to-leading order (NLO) inclusive multigluon 

spectrum was shown to take the form 



O = 



+ ^ E/ J2^^ J t«-fcAa • T„] [a+kXa ■ T,] I ,(2) 

A, a ^ J 

up to terms that do not contribute at leading logarithmic accuracy. The fields 
and 

^±k\ai-^) small quantum fluctuations propagating over the classi- 
cal field A'^. The former has a vanishing boundary condition in the remote past 
and its evolution is driven by a source term consisting of a 1-loop tadpole graph. 
The latter evolves without any source term, but its boundary condition in the 
remote past is a plane wave T"e^ ^±tk-x ^ ^ denote the polarization, 
color and momentum of the initial fluctuation). Here, d^u is the measure on S. 



^With this convention, y is the rapidity measured from the beam rapidity and y = 
corresponds to including only valence partons that have < < 1/P^. 



4 



Further, O^^ is a functional of the value of the classical field A'^ on the surface 
E and the differential operator acting on O^o generator of shifts of 

the value of A'^ on S. The only part of the expression for this operator that is 
important for computing leading logarithmic contributions i^ 

a.T.^a,(f^-X) ,,^(^i^.(,)) ^ (3) 

with rii^2 defined as in eq. (P) with y = ln(P±/A±). 

The next step in evaluating the NLO corrections is to integrate over the 
momentum k in eq. ([2]). One integral appears explicitly in the term that has 
two operators T, and another momentum integral is hidden in the source term of 
the fluctuation field . Since in the CGC effective theory the modes described 
as fields are bounded by < A* , these longitudinal momentum integrals have 
an upper bound. We shall compute only the contribution of modes in the small 
slices A'"*" < fc+ < A+ and A'~ < k~ < A~. By integrating out the field modes 
in these slices, one is going from the original CGC EFT to a new CGC EFT. 
The latter differs from the former because it has an additional layer of (slower) 
sources while its field modes now extend only up to smaller cutoffs A'='= . 

Using the results of paper I, we have to leading logarithmic accuracy. 



A, a 



E 



In (ttt) + In (^) Ha- . (4) 

A'±</c±<A± / / 

In this equation, Ha± are the JIMWLK Hamiltonians of the right and left 
moving nuclei respectively, at the scales A^ . For the nucleus moving in the +z 
direction, the explicit form of the JIMWLK Hamiltonian is 

2 J dAt{e-,y^) SAj{e-,xj_) 

with = 1/A+ and wher^ 

, N i /■ (a;l--»i)(j/l-txi) 
ryi (x^,y^)-^J ^^^^^ (a,^ - _ „^)2 

J ab 

There is a similar definition for the second nucleus moving in the — z direction. 



■^The complete operator a ■ T„ is made of three terms, but two of thom (not written 
explicitly here) do not provide any leading logarithmic contributions. 

Wilson lines without a rapidity argument are defined as in eq. {ij with y = ln(P±/A±) 
— they integrate over all the sources of the CGC EFT down to the cutoff A± . Moreover, the 
derivatives in the JIMWLK Hamiltonian are with respect to the slowest sources of the EFT. 
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An important point to note here is that the relation in eq. ([3]) is a property 
of the operator enclosed in the curly brackets of eq. ([2]), regardless of the details 
of the observable O under consideration. The second remarkable property of 
this result is that the leading logarithms can be assigned to one or the other of 
the two nuclei, without any mixing that would violate factorization. 

The expectation value of O, at NLO, can be represented in the CGC effective 
theory as 

0.o+(^N.o)= / [Dn^{y,x^)Dn2{y,x^)\ 

xW^^p^{y,x^)\W^_p2{y,x^)\ [O^o + O^^^], (7) 

where W^^ [r2i,2(y, x±)\ are the functional probability distributions for the Wil- 
son line configurations r2i,2(y, 

Inserting eq. ^ in eq. ^ and substituting the resulting expression on the 
r.h.s. of eq. ([7]), one can perform an integration by parts such that 7iA± 
operates on the distributions W^^ . Let us denotcj^ 

W^,^[n[{y,x^)] EE [l+ln(^|J^HA+]l^^Jf]i(y,aJx)] • (8) 

dy 

This equation is the infinitesimal form of the JIMWLK evolution equation, 
where n'l is the Wilson line corresponding to a source distribution p'l = pi+ 5pi 
and 5pi has support only in the interval [y, y + dy\. 

The argument of W^,^ in eq. ([8]) extends one step further in rapidity than 
the argument of - it is defined over the range < y < ln(_P+/A'+) and is 
hence a Wilson line in the CGC EFT. Thus our derivation proves that 

)-{o.o)' ■ (9) 

A'±<fe±<A± 

The prime on the r.h.s. indicates that the average is performed with a probabil- 
ity distribution corresponding to the Wilson lines of the CGC effective theory. 
In other words, this identity states that the classical expectation value of O 
in the original EFT, corrected by quantum fluctuations in a small slice of field 
modes, can be expressed as a purely classical expectation value in a new EFT 
with a lower cutoff and with a distribution of Wilson lines evolved according to 
eq. dH). 

Equation ([9]) describes how to resum the leading logarithmic quantum correc- 
tions in a small slice of longitudinal momentum. Successive leading logarithmic 
contributions down to fc^ = are obtained by repeating this elementary step 
infinitely many times while letting the thickness of the slices go to zero. One 



5 The JIMWLK Hamiltonian is Hermitian. 

''An identical analysis also applies to the second nucleus. 
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then obtains 



lo) = / [D^i{y,x^)DQ2{y,x^)\ Wp^{y,x^)\Wp2{v,x^)\ O^^,{l0) 

\ ' LLog J 

where W = hmA±^ol^A±" This expression is the central resuh of this paper. 
Eq. PU)) shows that aU the leading logarithms of rapidity, whether they cor- 
respond to the rapidity intervals between the nuclei and the tagged gluons or 
between the various tagged gluons, can be absorbed into the probability dis- 
tributions W for the trajectories of Wilson lines of the two projectiles. This 
formula applies to any inclusive observable for which eq. (0) is valid, regardless 
of whether the observable is local in rapidity or not. Because our result contains 
an average over y-dependent "trajectories" of Wilson lines, rather than an av- 
erage over Wilson lines at a given fixed rapidity, it contains a lot of information 
about multigluon correlations at different rapidities. 



3 One- and two-gluon inclusive spectra 




Figure 2: Diagrammatic representation of the various building blocks in the 
factorized formula for the inclusive single gluon spectrum. The lower part of 
the figure, representing nucleus 2, is made up of identical building blocks. 

We will now extract from our general result in eq. (jlOp . expressions for 
single and double inclusive gluon spectra. The single inclusive gluon spectrum 
dN\/d?p at LO depends only on the Wilson lines ^i.2{y,x±) at the rapidity 
y = Up of the produced gluon and not on the whole rapidity range as in eq. . 
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Therefore, we can simplify eq. ([TU| by inserting the identity 

1 = J [DUi,2ix^)] s[Ui,2ixi_)^niMyp,x±)] (11) 

and by defining the corresponding probabihty distributions for configurations 
of Wilson lines at the rapidity j/p 

Zy^[Ui,2ix^)] = J [Dni^2{y,x^)] w[ni^2iy,x^)] s[Ui^2{x^)-niMyp^x^)] . 

(12) 

One then obtains the all order leading log result for the single inclusive gluon 
spectrum in the following form 



/ dA^i 
\(Pp^dyp 



[DUiix^)DU2ix^)] Zy^ [Ui] Zy^ [U2] 



dN,[Ui,U2] 



d'^p^dyp 



(13) 



Note that the distribution Zy^[U] obeys the JIMWLK equation, 

dyMu]^ny^Zy^[u], (14) 

which must be supplemented by an initial condition at a rapidity close to the 
fragmentation region of the projectiles. Eq. is illustrated in the figure O 

At Leading Order, the inclusive two-gluon spectrum is simply the discon- 
nected product of two single gluon spectra |14I17| each of which depends on Wil- 
son lines at the rapidity of the corresponding gluon. Using eq. pT|) in eq. ifTO)) 
(now we need to insert four such delta functions) , one obtains for the resummed 
two-gluon spectrum the expression 



dA^2 



d^Pj_dypd'^q^dyq 



= J [DUf{x^)DUPix^)DUf{x^)DU^{x^)] 



dm[uiu^] 



- Zyp,yg [t^i ' ^i] Zyp,v, [^2 > U^] ^2p^^y^ 



dN,[UlU!i 



d^q^dy. 



,(15) 



where we have introduced the double probability distribution of Wilson lines as 



7 IttP Til 



'Dni^2iy,x^)] W[ni^2iy,x^)] 



x5[Ul2{x±)~fll.2{yp,X±)]S[Ul2{x±)~ni^2{yq,X±)] 



(16) 



This double distribution obeys the JIMWLK equation with respect to the largest 
of the two rapidities, 

if y, > yp, dy^Zy^.y^ [[/^ U ^] = Uy, Zy^.y^ [ [/^ , f/^] , (l7) 

with the boundary condition 

lim Zy^.y^ [U^, U^] = Zy^ [C/P] 5[U^ - [/'] . (18) 
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Alternately, this double distribution can be expressed as 



Zy.,y. [t^^ = {U\ U^^\ Z,^ [[/^'] , (19) 

where the Green's function Gy^_y^ [t/'', C/^] satisfies the JIMWLK equatiorQ 

dyfiy,,y, [U\ UP] = Hy, Gy^,y^ [U\ U^] , (20) 

with the initial condition 



Vq^Vp 



(21) 



This Green's function describes multigluon evolution, between two tagged glu- 
ons, in the presence of strong color sources from the projectiles. 

Our result for the double inclusive gluon spectrum, to leading logarithmic 
accuracy, can thus be expressed as follows 



diV, 



d2pj_dypd2qt^dy, 



= j [DUl{x^)DUl^{x^)DUl{x^)DUl{x^)\ 



'-Zy, [Ul] Gy^,y^ [C/f , Ul] Zy^ [U^] Gy^^y^ [[/|, U^] 



dN^[UlUl] 



d^Pxdyp 



d7Vi[[/f,C/|] 



(22) 



Eq. dm generahzes the result in paper II ~ that result, as implied by eq. (PT|) . 
is recovered when the rapidities of the two gluons are close to each other. Our 
formula for the two-gluon spectrum in eq. (|22|) is illustrated in fig. ([3]). By using 
eq. (jlOp . it is straightforward to write down similar formulae for higher gluon 
correlations. 

Factorization is obviously manifest in eq. (jlOp . It is slightly less apparent in 
eq. (|22p because there are more factors in the integrand. However, both the Z- 
functionals and the G Green's functionals are universal objects that describe the 
partonic content of a nucleus at high energy. The part of the integrand which is 
specific to the observable under consideration is relegated to the factors dNi/(Pp 
and dNi/tfq. 

Note also that these distributions are far more general than the fc_L -de- 
pendent correlators of Wilson lines that are often discussed in the literature 
of high energy QCD. The latter appear in a form of factorization which goes 
under the rubric of fcx-factorization. This type of factorization is formulated 
in terms of single gluon distributions, but is known to be violated |16| for 2- 
gluon correlations in collisions involving at least one saturated projectile. It is 
in this context that one should interpret the results of [TS] which concludes that 
factorization is violated. In stark contrast, our universal density functionals and 



'^We use here, and previously, the fact that the JIMWLK Hamiltonian acts only on objects 
at equal or higher rapidities. 
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from Yq to Yp 




Figure 3: Diagrammatic representation of the various building blocks in the 
factorized formula for the inclusive 2-gluon spectrum. As in the previous figure, 
the corresponding evolution from nucleus 2 at the bottom of the figure is not 
shown explicitly. 



Green's functions contain all the relevant information on rapidity dependent n- 
gluon correlations. In fact, our factorization result is a general consequence of 
causality and for this reason should even be valid beyond leading log accuracy. 

In practice, solving the JIMWLK equation to compute the r.h.s of eq. ([^ 
is more conveniently achieved by writing this equation as a Langevin equa- 
tion for Wilson lines living on the SU(3) group manifold [TB]. This stochas- 
tic formulation was implemented in the only extant numerical study of the 
JIMWLK equation [19] . Because solving the JIMWLK equation can be numer- 
ically challenging, a simpler formulation of high energy evolution is provided by 
the Balitsky-Kovchegov (BK) equation [20], which is a closed form mean field 
simplification of the JIMWLK expression for 2-point Wilson line correlators in 
the CGC E|. The BK equation corresponds to a nonlocal Gaussian form of the 
Zy^ , Zy^ functionals [21] (when expressed in terms of the color source distri- 
butions) with a variance fJ,\{yp^q, x±). Because the weight functionals at both 
Hp and Hq have this Gaussian form, the Green's function Gy^^y^ connecting the 
two must also be a Gaussian, whose variance can be determined from numerical 
solutions of the BK equation. Therefore, quantitative results for eq. ([^ can 
be obtained within this BK framework. These will be discussed in future work. 

*The BK equation is valid in the large A^c limit for large nuclei. 
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4 Dilute-dense limit 



The results obtained here are vahd for the coUision of two dense projectiles whose 
color charge densities of both are given by pi^2 ~ ■ This is the case in large 
nuclei or in nucleons at very high energies. An interesting question is whether 
decreasing the magnitude of pi and/or p2 in the formulas we have obtained so 
far gives the correct answer for dilute-dense or dilute-dilute collisions. Before 
going further, let us note that the permitted range for dilute to dense color 
charge densities is between g and g~^. The upper value, assumed in our study, 
corresponds to a fully saturated projectile. The lower value corresponds to a 
very dilute projectile whose parton density is of order unity ensuring that its 
density of color charge is proportional to g. 

The answer to the question posed is affirmative for the single inclusive gluon 
spectrum (eq. (jl3[) V By taking the limit p2 ^ g in this formula, one recov- 
ers immediately the well known result for the single inclusive spectrum in pA 
collisions and likewise for pp collisions when we let both pi^2 become of order g. 




Figure 4: Order of magnitude of various contributions to the 2-gluon spectrum as 
a function of the color charge density p2 in the small projectile (the color charge 
density in the large projectile is held fixed pi = 0{g~^)). Large logarithms of 
the energy, which become relevant in the leading logarithmic resummation, are 
not considered here. 

However, taking this limit in the 2-gluon spectrum (eq. ((22|) ) does not lead to 
the correct results for the inclusive 2-gluon spectrum in pA or pp collisions [16] . 
The reason of this discrepancy is that the corresponding power counting for the 
two gluon spectrum is very different for "dense" color sources ^ g^^ relative 
to the case of "dilute" sources ~ 5. In the power counting for dense sources 
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~ g~^, certain graphs are suppressed that would also be leading graphs for 
dilute sources ~ g. This is illustrated in figure IH where we have displayed the 
order of magnitude of three different contributions to the 2-gluon spectrum, as 
a function of p2- (For pA coUisions, pi is fixed to be of order g~^.) One sees 
that in the dense case, only one of these graphs is important, while they are 
all important in the dilute limit. Since eq. (|22p is obtained by assuming dense 
projectiles, it contains only the disconnected graph and misses the other two 

- this implies that eq. is not the complete answer in the dilute limit. As 
noted, this subtlety affects correlations between two or more gluons but not the 
single gluon spectrum. 

The previous discussion only deals with contributions at Leading Order. 
However, a similar discrepancy between the power counting in the dense-dense 
and dilute-dense limits occurs in the resummation of the leading logarithmic 
contributions. There, one sees that the operators included in the JIMWLK 
Hamiltonian are not the only ones that contribute at leading log accuracy in the 
dilute regime. For instance, when p ~ 5, an operator of the form g'^fP{S/S'p)'^ has 
the same order of magnitude as the operators in the JIMWLK Hamiltonian in 
the dilute limit (but is totally suppressed in the dense regime). Such an operator, 
with a prefactor of order and four derivatives with respect to the color source, 
corresponds to "pomeron splittings" in the energy evolution [16l22l23l24l25j 

- while the JIMWLK evolution only has "pomeron mergings" (because the 
number of p's in the prefactor is always equal to or greater than the number of 
derivatives). In principle, one would like a formalism where both limits contain 
the right physics. Unlike previous works which address the full S-matrix for 
high energy scattering, our focus will be on the more limited goal of computing 
inclusive gluon spectra in dense-dilute collisions. We believe that substantial 
progress in this direction is feasible and will further address this topic in a 
future publication. 

5 Conclusions 

In this paper, we obtained in eq. ([TO)) a general result for inclusive n-gluon 
production at arbitrary rapidities in the collision of two dense projectiles (such as 
heavy nuclei) with charge densities given by pi.2 ^ g^^- The result is expressed 
in terms of universal M^-density matrix functionals which contain information on 
n-gluon correlations in the wavefunctions of the dense projectiles. Our formalism 
is strictly valid in the leading logarithmic approximation in x. We anticipate, 
on the basis of simple causality arguments, that the structure of our result will 
hold beyond leading logarithmic accuracy. 

We explicitly wrote down the corresponding expressions for single and double 
inclusive gluon production with arbitrarily large rapidity separations between 
tagged gluons. Until this point, there was no microscopic QCD based formalism 
that allows the computation of the near side ridge correlations in nucleus-nucleus 
collisions when the rapidity separation between the measured particles is of the 
order ofl/as or more. Our formalism fills this gap and allows for future quanti- 
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tative comparisons and predictions for the rapidity dependence of the the ridge 
hke structures observed in central nucleus-nucleus collisions at RHIC and in 
future at the LHC. At the LHC, one may have the possibility of studying such 
structures that may span 6-10 units in rapidity. Such long range correlations 
therefore open a new window on the study of multiparton correlations in QCD 
as well as the provide a "chronometer" of the strong field initial "Glasma" stage 
of heavy ion collisions. 

When both projectiles are dense, we argued that only "Pomeron mergings" 
that are fully included in the JIMWLK Hamiltonian are relevant for inclu- 
sive gluon production. When one or both of the projectiles becomes dilute, 
our power counting suggests that "Pomeron splitting" contributions, become 
equally important as the "merging" contributions for correlations involving two 
or more gluons. Because these are not included in the JIMWLK formalism, they 
cannot be obtained by taking a naive low density limit of the dense-dense for- 
malism discussed in this paper. A smooth interpolation from the dilute-dense 
to the dense-dense limits for multigluon inclusive distributions requires that 
we first compute corrections to the JIMWLK Hamiltonian in the dilute-dense 
limit. While there have been several such studies in the context of the S-matrix 
for high energy scattering, they are in their infancy for inclusive multigluon 
production These studies will be important for extending our studies for 
nucleus-nucleus collisions to asymmetric systems such as high energy proton- 
nucleus collisions. This work is in progress and will be addressed in a future 
publication. 

Acknowledgements 

We thank the Center for Theoretical Sciences of the Tata Institute for Fun- 
damental Research for their support during the program "Initial Conditions 
in Heavy Ion Collisions". R.V.'s research is supported by the US Department 
of Energy under DOE Contract No. DE-AC02-98CH10886. F.G.'s work is 
supported in part by Agence Nationale de la Recherche via the programme 
ANR-06-BLAN-0285-01. 

References 

[1] J. Putschke, J. Phys. G34, S679 (2007). 



[2] STAR, M. Daugherity, |arXiv:0806.2121 [nucl-ex 



[3] PHOBOS, B. Alver et al, J. Phys. G35, 104080 (2008). 

[4] STAR, B.K. Srivastava, Int. J. Mod. Phys. E16, 3371 (2007). 



5] E. lancu and R. Venugopalan, arXiv:hep-ph/0303204[ H. Weigert, Prog. 



Part. Nucl. Phys. 55, 461 (2005). 



13 



[6] L. D. McLerran and R. Venugopalan, Phys. Rev. D49, 2233 (1994); ibid. 
3352 (1994); ibtd. D50, 2225 (1994). 

[7] J. Jalilian-Marian, A. Kovner, L. D. McLerran and H. Weigert, Phys. 
Rev. D55, 5414 (1997); J. Jalilian-Marian, A. Kovner, A. Leonidov and 
H. Weigert, Nucl. Phys. B504, 415 (1997); ibid. Phys. Rev. D59, 014014 
(1998); E. lancu, A. Leonidov and L. D. McLerran, Nucl. Phys. A692, 
583 (2001); ibid. Phys. Lett. B510, 133 (2001); E. Ferreiro, E. lancu, 
A. Leonidov and L. McLerran, Nucl. Phys. A703, 489 (2002). 

T. Lappi and L. McLerran, Nucl. Phys. A772, 200 (2006). 

A. Kovner, L. D. McLerran and H. Weigert, Phys. Rev. D52, 3809 (1995). 

D. Kharzeev, A. Krasnitz and R. Venugopalan, Phys. Lett. B545, 298 
(2002). 

A. Dumitru, F. Gelis, L. McLerran and R. Venugopalan, Nucl. Phys. A810, 
91 (2008); 



S. Gavin, L. McLerran and G. MoscheUi, |arXiv:0806.4718 [nucl-th f, 

F. Gelis, T. Lappi and R. Venugopalan, Phys. Rev. D78, 054019 (2008). 

F. Gelis, T. Lappi and R. Venugopalan, Phys. Rev. D78, 054020 (2008). 

M. Ciafaloni and G. Marchesini, Nucl. Phys. B109, 261 (1976); M.A. 
Braun, Eur. Phys. J. C48, 501 (2006); ibid. arXiv:0801.0493. 

J. Jalihan-Marian and Y. V. Kovchegov, Phys. Rev. D70, 114017 (2004). 

D. Kharzeev, E. Levin and L. McLerran, Nucl. Phys. A748, 627 (2005); 
N. Armesto, L. McLerran and C. Pajares, Nucl. Phys. A781, 201 (2007). 

J.-P. Blaizot, E. lancu and H. Weigert, Nucl. Phys. A713, 441 (2003). 

K. Rummukaincn and H. Weigert, Nucl. Phys. A739, 183 (2004). 

L Bahtsky, Nucl. Phys. B463, 99 (1996); Y. V. Kovchegov, Phys. Rev. 
D60, 034008 (1999). 

H. Fujii, F. Gehs and R. Venugopalan, Nucl. Phys. A780, 146 (2006). 

A. H. Mueller, A. I. Shoshi and S. M. H. Wong, Nucl. Phys. B 715 (2005) 
440. 

E. lancu and D. N. Triantafyllopoulos, Nucl. Phys. A 756 (2005) 419. 
A. Kovner and M. Lubhnsky, Phys. Rev. D 71 (2005) 085004. 

A. Kovner and M. Lubhnsky, Phys. Rev. D 72 (2005) 074023. 
A. Kovner and M. Lubhnsky, JHEP 11, 083 (2006). 



14 



